Unveiling the link between physical activity levels and dementia risk: Insights from the UK Biobank study.

BACKGROUND: There is limited information on the mixture effect and weights of light physical activity (LPA), moderate physical activity (MPA), and vigorous physical activity (VPA) on dementia risk. METHODS: A prospective cohort study was conducted based on the UK Biobank dataset. We included participants aged at least 45 years old without dementia at baseline between 2006-2010. The weighted quantile sum regression was used to explore the mixture effect and weights of three types of physical activity on dementia risk. RESULTS: This study includes 354,123 participants, with a mean baseline age of 58.0-year-old and 52.4 % of female participants. During a median follow-up time of 12.5 years, 5,136 cases of dementia were observed. The mixture effect of LPA, MPA, and VPA on dementia was statistically significant (ß: -0.0924, 95 % Confidence Interval (CI): (-0.1402, -0.0446), P < 0.001), with VPA (weight: 0.7922) contributing most to a lower dementia risk, followed by MPA (0.1939). For Alzheimer's disease, MPA contributed the most (0.8555); for vascular dementia, VPA contributed the most (0.6271). CONCLUSION: For Alzheimer's disease, MPA was identified as the most influential factor, while VPA stood out as the most impactful for vascular dementia.

Somatomotor-visual resting state functional connectivity increases after 2 years in the UK Biobank longitudinal cohort.

Purpose: Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, in which high connectivity among all brain regions changes to a more modular structure with maturation. We examine FC changes in older adults after 2 years of aging in the UK Biobank (UKB) longitudinal cohort. Approach: We process fMRI connectivity data using the Power264 atlas and then test whether the average internetwork FC changes in the 2722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, independent component analysis (ICA)-based FC to determine which of a longitudinal scan pair is older. Finally, we investigate cross-sectional FC changes as well as differences due to differing scanner tasks in the UKB, Philadelphia Neurodevelopmental Cohort, and Alzheimer's Disease Neuroimaging Initiative datasets. Results: We find a 6.8% average increase in somatomotor network (SMT)-visual network (VIS) connectivity from younger to older scans (corrected p<10-15) that occurs in male, female, older subject (>65 years old), and younger subject (<55 years old) groups. Among all internetwork connections, the average SMT-VIS connectivity is the best predictor of relative scan age. Using the full FC and a training set of 2000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC. Conclusions: We conclude that SMT-VIS connectivity increases with age in the UKB longitudinal cohort and that resting state FC increases with age in the UKB cross-sectional cohort.

Multi-View Variational Autoencoder for Missing Value Imputation in Untargeted Metabolomics.

Background: Missing data is a common challenge in mass spectrometry-based metabolomics, which can lead to biased and incomplete analyses. The integration of whole-genome sequencing (WGS) data with metabolomics data has emerged as a promising approach to enhance the accuracy of data imputation in metabolomics studies. Method: In this study, we propose a novel method that leverages the information from WGS data and reference metabolites to impute unknown metabolites. Our approach utilizes a multi-view variational autoencoder to jointly model the burden score, polygenetic risk score (PGS), and linkage disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature extraction and missing metabolomics data imputation. By learning the latent representations of both omics data, our method can effectively impute missing metabolomics values based on genomic information. Results: We evaluate the performance of our method on empirical metabolomics datasets with missing values and demonstrate its superiority compared to conventional imputation techniques. Using 35 template metabolites derived burden scores, PGS and LD-pruned SNPs, the proposed methods achieved R2-scores > 0.01 for 71.55% of metabolites. Conclusion: The integration of WGS data in metabolomics imputation not only improves data completeness but also enhances downstream analyses, paving the way for more comprehensive and accurate investigations of metabolic pathways and disease associations. Our findings offer valuable insights into the potential benefits of utilizing WGS data for metabolomics data imputation and underscore the importance of leveraging multi-modal data integration in precision medicine research.

Multi-view information fusion using multi-view variational autoencoder to predict proximal femoral fracture load.

Background: Hip fracture occurs when an applied force exceeds the force that the proximal femur can support (the fracture load or "strength") and can have devastating consequences with poor functional outcomes. Proximal femoral strengths for specific loading conditions can be computed by subject-specific finite element analysis (FEA) using quantitative computerized tomography (QCT) images. However, the radiation and availability of QCT limit its clinical usability. Alternative low-dose and widely available measurements, such as dual energy X-ray absorptiometry (DXA) and genetic factors, would be preferable for bone strength assessment. The aim of this paper is to design a deep learning-based model to predict proximal femoral strength using multi-view information fusion. Results: We developed new models using multi-view variational autoencoder (MVAE) for feature representation learning and a product of expert (PoE) model for multi-view information fusion. We applied the proposed models to an in-house Louisiana Osteoporosis Study (LOS) cohort with 931 male subjects, including 345 African Americans and 586 Caucasians. We performed genome-wide association studies (GWAS) to select 256 genetic variants with the lowest p-values for each proximal femoral strength and integrated whole genome sequence (WGS) features and DXA-derived imaging features to predict proximal femoral strength. The best prediction model for fall fracture load was acquired by integrating WGS features and DXA-derived imaging features. The designed models achieved the mean absolute percentage error of 18.04%, 6.84% and 7.95% for predicting proximal femoral fracture loads using linear models of fall loading, nonlinear models of fall loading, and nonlinear models of stance loading, respectively. Conclusion: The proposed models are capable of predicting proximal femoral strength using WGS features and DXA-derived imaging features. Though this tool is not a substitute for predicting FEA using QCT images, it would make improved assessment of hip fracture risk more widely available while avoiding the increased radiation exposure from QCT.

Metabolomic profiles associated with physical activity in White and African American adult men.

BACKGROUND: Physical activity (PA) is associated with various health benefits, especially in improving chronic health conditions. However, the metabolic changes in host metabolism in response to PA remain unclear, especially in racially/ethnically diverse populations. OBJECTIVE: This study is to assess the metabolic profiles associated with the frequency of PA in White and African American (AA) men. METHODS: Using the untargeted metabolomics data collected from 698 White and AA participants (mean age: 38.0±8.0, age range: 20-50) from the Louisiana Osteoporosis Study (LOS), we conducted linear regression models to examine metabolites that are associated with PA levels (assessed by self-reported regular exercise frequency levels: 0, 1-2, and ≥3 times per week) in White and AA men, respectively, as well as in the pooled sample. Covariates considered for statistical adjustments included race (only for the pooled sample), age, BMI, waist circumstance, smoking status, and alcohol drinking. RESULTS: Of the 1133 untargeted compounds, we identified 7 metabolites associated with PA levels in the pooled sample after covariate adjustment with a false discovery rate of 0.15. Specifically, compared to participants who did not exercise, those who exercised at a frequency ≥3 times/week showed higher abundances in uracil, orotate, 1-(1-enyl-palmitoyl)-2-oleoyl-GPE (P-16:0/18:1) (GPE), threonate, and glycerate, but lower abundances in salicyluric glucuronide and adenine in the pooled sample. However, in Whites, salicyluric glucuronide and orotate were not significant. Adenine, GPE, and threonate were not significant in AAs. In addition, the seven metabolites were not significantly different between participants who exercised ≥3 times/week and 1-2 times/week, nor significantly different between participants with 1-2 times/week and 0/week in the pooled sample and respective White and AA groups. CONCLUSIONS: Metabolite responses to PA are dose sensitive and may differ between White and AA populations. The identified metabolites may help advance our knowledge of guiding precision PA interventions. Studies with rigorous study designs are warranted to elucidate the relationship between PA and metabolites.

Gut microbiota impacts bone via Bacteroides vulgatus-valeric acid-related pathways.

Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

omicsMIC: a Comprehensive Benchmarking Platform for Robust Comparison of Imputation Methods in Mass Spectrometry-based Omics Data.

Mass spectrometry is a powerful and widely used tool for generating proteomics, lipidomics, and metabolomics profiles, which is pivotal for elucidating biological processes and identifying biomarkers. However, missing values in spectrometry-based omics data may pose a critical challenge for the comprehensive identification of biomarkers and elucidation of the biological processes underlying human complex disorders. To alleviate this issue, various imputation methods for mass spectrometry-based omics data have been developed. However, a comprehensive and systematic comparison of these imputation methods is still lacking, and researchers are frequently confronted with a multitude of options without a clear rationale for method selection. To address this pressing need, we developed omicsMIC (mass spectrometry-based omics with Missing values Imputation methods Comparison platform), an interactive platform that provides researchers with a versatile framework to simulate and evaluate the performance of 28 diverse imputation methods. omicsMIC offers a nuanced perspective, acknowledging the inherent heterogeneity in biological data and the unique attributes of each dataset. Our platform empowers researchers to make data-driven decisions in imputation method selection based on real-time visualizations of the outcomes associated with different imputation strategies. The comprehensive benchmarking and versatility of omicsMIC make it a valuable tool for the scientific community engaged in mass spectrometry-based omics research. OmicsMIC is freely available at https://github.com/WQLin8/omicsMIC.

Identifiability in Functional Connectivity May Unintentionally Inflate Prediction Results.

Functional magnetic resonance (fMRI) is an invaluable tool in studying cognitive processes in vivo. Many recent studies use functional connectivity (FC), partial correlation connectivity (PC), or fMRI-derived brain networks to predict phenotypes with results that sometimes cannot be replicated. At the same time, FC can be used to identify the same subject from different scans with great accuracy. In this paper, we show a method by which one can unknowingly inflate classification results from 61% accuracy to 86% accuracy by treating longitudinal or contemporaneous scans of the same subject as independent data points. Using the UK Biobank dataset, we find one can achieve the same level of variance explained with 50 training subjects by exploiting identifiability as with 10,000 training subjects without double-dipping. We replicate this effect in four different datasets: the UK Biobank (UKB), the Philadelphia Neurodevelopmental Cohort (PNC), the Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP), and an OpenNeuro Fibromyalgia dataset (Fibro). The unintentional improvement ranges between 7% and 25% in the four datasets. Additionally, we find that by using dynamic functional connectivity (dFC), one can apply this method even when one is limited to a single scan per subject. One major problem is that features such as ROIs or connectivities that are reported alongside inflated results may confuse future work. This article hopes to shed light on how even minor pipeline anomalies may lead to unexpectedly superb results.

Somatomotor-Visual Resting State Functional Connectivity Increases After Two Years in the UK Biobank Longitudinal Cohort.

Purpose: Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, where high connectivity among all brain regions changes to a more modular structure with maturation. In this work, we examine changes in FC in older adults after two years of aging in the UK Biobank longitudinal cohort. Approach: We process data using the Power264 atlas, then test whether FC changes in the 2,722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, ICA-based FC to determine which of a longitudinal scan pair is older. Results: We find a 6.8% average increase in SMT-VIS connectivity from younger to older scan (from ρ = 0.39 to ρ = 0.42) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all inter-network connections, this average SMT-VIS connectivity is the best predictor of relative scan age, accurately predicting which scan is older 57% of the time. Using the full FC and a training set of 2,000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC. Conclusions: We conclude that SMT-VIS connectivity increases in the longitudinal cohort, while resting state FC increases generally with age in the cross-sectional cohort. However, we consider the possibility of a change in resting state scanner task between UKB longitudinal data acquisitions.

Somatomotor-Visual Resting State Functional Connectivity Increases After Two Years in the UK Biobank Longitudinal Cohort.

Purpose: Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, where high connectivity among all brain regions changes to a more modular structure with maturation. In this work, we examine changes in FC in older adults after two years of aging in the UK Biobank longitudinal cohort. Approach: We process data using the Power264 atlas, then test whether FC changes in the 2,722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, ICA-based FC to determine which of a longitudinal scan pair is older. Results: We find a 6.8% average increase in SMT-VIS connectivity from younger to older scan (from ρ=0.39 to ρ=0.42) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all inter-network connections, this average SMT-VIS connectivity is the best predictor of relative scan age, accurately predicting which scan is older 57% of the time. Using the full FC and a training set of 2,000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC. Conclusions: We conclude that SMT-VIS connectivity increases in the longitudinal cohort, while resting state FC increases generally with age in the cross-sectional cohort. However, we consider the possibility of a change in resting state scanner task between UKB longitudinal data acquisitions.

Identification of potential genetic causal variants for obesity-related traits using statistical fine mapping.

BACKGROUND: Obesity is highly influenced by heritability and variant effects. While previous genome-wide association studies (GWASs) have successfully identified numerous genetic loci associated with obesity-related traits [body mass index (BMI) and waist-to-hip ratio (WHR)], most causal variants remain unidentified. The high degree of linkage disequilibrium (LD) throughout the genome makes it extremely difficult to distinguish the GWAS-associated SNPs that exert a true biological effect. OBJECTIVE: This study was to identify the potential causal variants having a biological effect on obesity-related traits. METHODS: We used Probabilistic Annotation INTegratOR, a Bayesian fine-mapping method, which incorporated genetic association data (GWAS summary statistics), LD structure, and functional annotations to calculate a posterior probability of causality for SNPs across all loci of interest. Moreover, we performed gene expression analysis using the available public transcriptomic data to validate the corresponding genes of the potential causal SNPs partially. RESULTS: We identified 96 SNPs for BMI and 43 SNPs for WHR with a high posterior probability of causality (> 99%), including 49 BMI SNPs and 24 WHR SNPs which did not reach genome-wide significance in the original GWAS. Finally, we partially validated some genes corresponding to the potential causal SNPs. CONCLUSION: Using a statistical fine-mapping approach, we identified a set of potential causal variants to be prioritized for future functional validation and also detected some novel trait-associated variants. These results provided novel insight into our understanding of the genetics of obesity and also demonstrated that fine mapping may improve upon the results identified by the original GWASs.

Systematic metabolomic studies identified adult adiposity biomarkers with acetylglycine associated with fat loss in vivo.

Obesity is associated with various adverse health outcomes. Body fat (BF) distribution is recognized as an important factor of negative health consequences of obesity. Although metabolomics studies, mainly focused on body mass index (BMI) and waist circumference, have explored the biological mechanisms involved in the development of obesity, these proxy composite measures are not accurate and cannot reflect BF distribution, and thus may hinder accurate assessment of metabolic alterations and differential risk of metabolic disorders among individuals presenting adiposity differently throughout the body. Thus, the exact relations between metabolites and BF remain to be elucidated. Here, we aim to examine the associations of metabolites and metabolic pathways with BF traits which reflect BF distribution. We performed systematic untargeted serum metabolite profiling and dual-energy X-ray absorptiometry (DXA) whole body fat scan for 517 Chinese women. We jointly analyzed DXA-derived four BF phenotypes to detect cross-phenotype metabolite associations and to prioritize important metabolomic factors. Topology-based pathway analysis was used to identify important BF-related biological processes. Finally, we explored the relationships of the identified BF-related candidate metabolites with BF traits in different sex and ethnicity through two independent cohorts. Acetylglycine, the top distinguished finding, was validated for its obesity resistance effect through in vivo studies of various diet-induced obese (DIO) mice. Eighteen metabolites and fourteen pathways were discovered to be associated with BF phenotypes. Six of the metabolites were validated in varying sex and ethnicity. The obesity-resistant effects of acetylglycine were observed to be highly robust and generalizable in both human and DIO mice. These findings demonstrate the importance of metabolites associated with BF distribution patterns and several biological pathways that may contribute to obesity and obesity-related disease etiology, prevention, and intervention. Acetylglycine is highlighted as a potential therapeutic candidate for preventing excessive adiposity in future studies.

An autoencoder-based deep learning method for genotype imputation.

Genotype imputation has a wide range of applications in genome-wide association study (GWAS), including increasing the statistical power of association tests, discovering trait-associated loci in meta-analyses, and prioritizing causal variants with fine-mapping. In recent years, deep learning (DL) based methods, such as sparse convolutional denoising autoencoder (SCDA), have been developed for genotype imputation. However, it remains a challenging task to optimize the learning process in DL-based methods to achieve high imputation accuracy. To address this challenge, we have developed a convolutional autoencoder (AE) model for genotype imputation and implemented a customized training loop by modifying the training process with a single batch loss rather than the average loss over batches. This modified AE imputation model was evaluated using a yeast dataset, the human leukocyte antigen (HLA) data from the 1,000 Genomes Project (1KGP), and our in-house genotype data from the Louisiana Osteoporosis Study (LOS). Our modified AE imputation model has achieved comparable or better performance than the existing SCDA model in terms of evaluation metrics such as the concordance rate (CR), the Hellinger score, the scaled Euclidean norm (SEN) score, and the imputation quality score (IQS) in all three datasets. Taking the imputation results from the HLA data as an example, the AE model achieved an average CR of 0.9468 and 0.9459, Hellinger score of 0.9765 and 0.9518, SEN score of 0.9977 and 0.9953, and IQS of 0.9515 and 0.9044 at missing ratios of 10% and 20%, respectively. As for the results of LOS data, it achieved an average CR of 0.9005, Hellinger score of 0.9384, SEN score of 0.9940, and IQS of 0.8681 at the missing ratio of 20%. In summary, our proposed method for genotype imputation has a great potential to increase the statistical power of GWAS and improve downstream post-GWAS analyses.

Integration of the Human Gut Microbiome and Serum Metabolome Reveals Novel Biological Factors Involved in the Regulation of Bone Mineral Density.

While the gut microbiome has been reported to play a role in bone metabolism, the individual species and underlying functional mechanisms have not yet been characterized. We conducted a systematic multi-omics analysis using paired metagenomic and untargeted serum metabolomic profiles from a large sample of 499 peri- and early post-menopausal women to identify the potential crosstalk between these biological factors which may be involved in the regulation of bone mineral density (BMD). Single omics association analyses identified 22 bacteria species and 17 serum metabolites for putative association with BMD. Among the identified bacteria, Bacteroidetes and Fusobacteria were negatively associated, while Firmicutes were positively associated. Several of the identified serum metabolites including 3-phenylpropanoic acid, mainly derived from dietary polyphenols, and glycolithocholic acid, a secondary bile acid, are metabolic byproducts of the microbiota. We further conducted a supervised integrative feature selection with respect to BMD and constructed the inter-omics partial correlation network. Although still requiring replication and validation in future studies, the findings from this exploratory analysis provide novel insights into the interrelationships between the gut microbiome and serum metabolome that may potentially play a role in skeletal remodeling processes.

Pathway-based metabolomics study of sarcopenia-related traits in two US cohorts.

We aimed to validate two metabolites, aspartic acid and glutamic acid, which were associated with sarcopenia-related traits, muscle mass and strength, in our previous untargeted metabolomics study and to identify novel metabolites from five metabolic pathways involving these two metabolites. We included a discovery cohort of 136 white women aged 20-40 years (used for the previous untargeted metabolomics analysis) and a validation cohort of 174 subjects aged ≥ 60 years, including men and women of white and black. A targeted LC-MS assay successfully detected 12 important metabolites from these pathways. Aspartic acid was associated with muscle mass and strength in the discovery cohort, but not in the validation cohort. However, glutamic acid was associated with these sarcopenia traits in both cohorts. Additionally, N-acetyl-L-aspartic acid and carnosine were the newly identified metabolites that were associated with muscle strength in the discovery and validation cohorts, respectively. We did not observe any significant sex and race differences in the associations of these metabolites with sarcopenia traits in the validation cohort. Our findings indicated that glutamic acid might be consistently associated with sarcopenia-related traits across age, sex, and race. They also suggested that age-specific metabolites and metabolic pathways might be involved in muscle regulation.

A multiethnic whole genome sequencing study to identify novel loci for bone mineral density.

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis.

A transcriptome-wide association study to detect novel genes for volumetric bone mineral density.

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm2) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm3) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation.

Enhanced Identification of Novel Potential Variants for Appendicular Lean Mass by Leveraging Pleiotropy With Bone Mineral Density.

Strong relationships have been found between appendicular lean mass (ALM) and bone mineral density (BMD). It may be due to a shared genetic basis, termed pleiotropy. By leveraging the pleiotropy with BMD, the aim of this study was to detect more potential genetic variants for ALM. Using the conditional false discovery rate (cFDR) methodology, a combined analysis of the summary statistics of two large independent genome wide association studies (GWAS) of ALM (n = 73,420) and BMD (n = 10,414) was conducted. Strong pleiotropic enrichment and 26 novel potential pleiotropic SNPs were found for ALM and BMD. We identified 156 SNPs for ALM (cFDR <0.05), of which 74 were replicates of previous GWASs and 82 were novel SNPs potentially-associated with ALM. Eleven genes annotated by 31 novel SNPs (13 pleiotropic and 18 ALM specific) were partially validated in a gene expression assay. Functional enrichment analysis indicated that genes corresponding to the novel potential SNPs were enriched in GO terms and/or KEGG pathways that played important roles in muscle development and/or BMD metabolism (adjP <0.05). In protein-protein interaction analysis, rich interactions were demonstrated among the proteins produced by the corresponding genes. In conclusion, the present study, as in other recent studies we have conducted, demonstrated superior efficiency and reliability of the cFDR methodology for enhanced detection of trait-associated genetic variants. Our findings shed novel insight into the genetic variability of ALM in addition to the shared genetic basis underlying ALM and BMD.

Identification and Functional Characterization of Metabolites for Bone Mass in Peri- and Postmenopausal Chinese Women.

CONTEXT: Although metabolic profiles appear to play an important role in menopausal bone loss, the functional mechanisms by which metabolites influence bone mineral density (BMD) during menopause are largely unknown. OBJECTIVE: We aimed to systematically identify metabolites associated with BMD variation and their potential functional mechanisms in peri- and postmenopausal women. DESIGN AND METHODS: We performed serum metabolomic profiling and whole-genome sequencing for 517 perimenopausal (16%) and early postmenopausal (84%) women aged 41 to 64 years in this cross-sectional study. Partial least squares regression and general linear regression analysis were applied to identify BMD-associated metabolites, and weighted gene co-expression network analysis was performed to construct co-functional metabolite modules. Furthermore, we performed Mendelian randomization analysis to identify causal relationships between BMD-associated metabolites and BMD variation. Finally, we explored the effects of a novel prominent BMD-associated metabolite on bone metabolism through both in vivo/in vitro experiments. RESULTS: Twenty metabolites and a co-functional metabolite module (consisting of fatty acids) were significantly associated with BMD variation. We found dodecanoic acid (DA), within the identified module causally decreased total hip BMD. Subsequently, the in vivo experiments might support that dietary supplementation with DA could promote bone loss, as well as increase the osteoblast and osteoclast numbers in normal/ovariectomized mice. Dodecanoic acid treatment differentially promoted osteoblast and osteoclast differentiation, especially for osteoclast differentiation at higher concentrations in vitro (eg,10, 100 µM). CONCLUSIONS: This study sheds light on metabolomic profiles associated with postmenopausal osteoporosis risk, highlighting the potential importance of fatty acids, as exemplified by DA, in regulating BMD.

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis.

CONTEXT: Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms. OBJECTIVE: Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis. DESIGN, SETTING, AND PARTICIPANTS: Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings. RESULTS: In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1. CONCLUSIONS: Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.